KSNP: a fast de Bruijn graph-based haplotyping tool approaching data-in time cost.

Long reads that cover more variants per read raise opportunities for accurate haplotype construction, whereas the genotype errors of single nucleotide polymorphisms pose great computational challenges for haplotyping tools. Here we introduce KSNP, an efficient haplotype construction tool based on the de Bruijn graph (DBG). KSNP leverages the ability of DBG in handling high-throughput erroneous reads to tackle the challenges. Compared to other notable tools in this field, KSNP achieves at least 5-fold speedup while producing comparable haplotype results. The time required for assembling human haplotypes is reduced to nearly the data-in time.

A chromosomal-scale genome assembly of modern cultivated hybrid sugarcane provides insights into origination and evolution.

Sugarcane is a vital crop with significant economic and industrial value. However, the cultivated sugarcane's ultra-complex genome still needs to be resolved due to its high ploidy and extensive recombination between the two subgenomes. Here, we generate a chromosomal-scale, haplotype-resolved genome assembly for a hybrid sugarcane cultivar ZZ1. This assembly contains 10.4 Gb genomic sequences and 68,509 annotated genes with defined alleles in two sub-genomes distributed in 99 original and 15 recombined chromosomes. RNA-seq data analysis shows that sugar accumulation-associated gene families have been primarily expanded from the ZZSO subgenome. However, genes responding to pokkah boeng disease susceptibility have been derived dominantly from the ZZSS subgenome. The region harboring the possible smut resistance genes has expanded significantly. Among them, the expansion of WAK and FLS2 families is proposed to have occurred during the breeding of ZZ1. Our findings provide insights into the complex genome of hybrid sugarcane cultivars and pave the way for future genomics and molecular breeding studies in sugarcane.

The effects of ANRIL polymorphisms on colorectal cancer, tumor stage, and tumor grade among Iranian population.

BACKGROUND: Colorectal cancer (CRC) is a type of neoplasm, developing in the colon or rectum. The exact etiology of CRC is not well known, but the role of genetic, epigenetic, and environmental factors are established in its pathogenesis. Therefore, the aim of this research was to explore the effects of ANRIL polymorphisms on the CRC and its clinical findings. METHODS AND RESULTS: The peripheral blood specimens were collected from 142 CRC patients and 225 controls referred to Milad Hospital, Tehran, Iran. PCR- RFLP method was used to analyze ANRIL rs1333040, rs10757274 rs4977574, and rs1333048 polymorphisms. The ANRIL rs1333040 polymorphism was related to a higher risk of CRC in the co-dominant, dominant, and log-additive models. ANRIL rs10757274, rs4977574, and rs1333048 polymorphisms showed no effect on CRC susceptibility. The CGAA and TGGA haplotypes of ANRIL rs1333040/ rs10757274/ rs4977574/rs1333048 polymorphisms were associated with the higher and the lower risk of CRC respectively. The rs1333040 polymorphism was associated with higher TNM stages (III and IV). The frequency of ANRIL rs10757274 polymorphism was lower in CRC patients over 50 years of age only in the dominant model. In addition, the rs10757274 was associated with well differentiation in CRC patients. CONCLUSION: The ANRIL rs1333040 polymorphism was associated with a higher risk of CRC and higher TNM stages. ANRIL rs10757274 polymorphism was associated with the well-differentiated tumor in CRC.

Haplotype-resolved assembly of auto-polyploid genomes via combining Hi-C and gametic data.

Haplotype-resolved genome assembly plays a crucial role in understanding allele-specific functions. However, obtaining haplotype-resolved assembly for auto-polyploid genomes remains challenging. Existing methods can be classified into reference-based phasing, assembly-based phasing, and gamete binning. Nevertheless, there is a lack of cost-effective and efficient methods for haplotyping auto-polyploid genomes. In this study, we propose a novel phasing algorithm called PolyGH, which combines Hi-C and gametic data. We conducted experiments on tetraploid potato cultivars and divided the method into three steps. Firstly, gametic data was utilized to bin non-collapsed contigs, followed by merging adjacent fragments of the same type within the same contig. Secondly, accurate Hi-C signals related to differential genomic regions were acquired using unique k-mers. Finally, collapsed fragments were assigned to haplotigs based on combined Hi-C and gametic signals. Comparing PolyGH with Hi-C-based and gametic data-based methods, we found that PolyGH exhibited superior performance in haplotyping auto-polyploid genomes when integrating both data types. This approach has the potential to enhance haplotype-resolved assembly for auto-polyploid genomes.

Nanopore sequencing with T2T-CHM13 for accurate detection and preventing the transmission of structural rearrangements in highly repetitive heterochromatin regions in human embryos.

BACKGROUND: Structural rearrangements in highly repetitive heterochromatin regions can result in miscarriage or foetal malformations; however, detecting and preventing the transmission of these rearrangements has been challenging. Recently, the completion of sequencing of the complete human genome (T2T-CHM13) has made it possible to accurately characterise structural rearrangements in these regions. We developed a method based on T2T-CHM13 and nanopore sequencing to detect and block structural rearrangements in highly repetitive heterochromatin sequences. METHODS: T2T-CHM13-based "Mapping Allele with Resolved Carrier Status" was performed for couples who carry structural rearrangements in heterochromatin regions. Using nanopore sequencing and the T2T-CHM13 reference genome, the precise breakpoints of inversions and translocations close to the centromere were detected and haplotypes were constructed using flanking single-nucleotide polymorphisms (SNPs). Haplotype linkage analysis was then performed by comparing consistent parental SNPs with embryonic SNPs to determine whether the embryos carried hereditary inversions or balanced translocations. Based on copy number variation and haplotype linkage analysis, we transplanted normal embryos, which were further verified by an amniotic fluid test. RESULTS: To validate this approach, we used nanopore sequencing of families with inversions and reciprocal translocations close to the centromere. Using the T2T-CHM13 reference genome, we accurately detected inversions and translocations in centromeres, constructed haplotypes and prevented the transmission of structural rearrangements in the offspring. CONCLUSIONS: This study represents the first successful application of T2T-CHM13 in human reproduction and provides a feasible protocol for detecting and preventing the transmission of structural rearrangements of heterochromatin in embryos.

Population genetic structure of Randall's threadfin bream Nemipterus randalli in Indian waters based on mitochondrial and nuclear gene sequences.

Nemipterus randalli, commonly known as Randall's threadfin bream, is a commercially important marine finfish. Understanding its genetic structure is critical to effective management and conservation efforts. Previous investigations on population structure in this species were limited by geographic coverage. In this study, we utilized the mitochondrial Cytochrome b gene and nuclear Ribosomal protein gene intron Rp S7 sequences to investigate the population genetic structure, demography and genetic diversity of N. randalli along Indian waters. Our results revealed high haplotype diversity but low nucleotide diversity. AMOVA revealed that the variation among the population was highly significant. Hierarchical AMOVA provided further evidence of significant genetic differentiation between the west and east coasts, which was corroborated by the Bayesian tree and the median-joining network diagram. The mtDNA sequences revealed significant genetic structure between populations based on fixation index analysis following the isolation-by-distance model. Furthermore, the neutrality test and mismatch analysis suggest that N. randalli populations may have experienced a population expansion. However, nuclear marker RpS7, showed a high level of polymorphism, which obscured the population structuring observed with the mitochondrial marker. Consequently, concordant results were not obtained when comparing the mitochondrial and nuclear DNA sequences. The strong genetic differentiation between the east and west coast observed using mitochondrial marker could be attributed to a combination of geographic and environmental factors. These findings lay the groundwork for developing effective conservation and management strategies for N. randalli, considering its genetic structure.

IGF2BP2 and IGFBP3 Genotypes, Haplotypes, and Genetic Models Studies in Polycystic Ovary Syndrome.

BACKGROUND: Insulin resistance has been correlated with the genetic diversity within the insulin-like binding proteins genes. Moreover, insulin resistance is one of the key characteristics of the widespread reproductive endocrine condition known as polycystic ovarian syndrome (PCOS). Hence, this study is aimed to determine the association between IGFBP3 and IGF2BP2 gene variants and PCOS risk. METHODS: A total of 300 subjects (150 PCOS cases diagnosed based on Rotterdam ESHRE/ASRM consensus criteria and 150 healthy subjects) were recruited in this case-control cross-sectional study. Tetra-primer amplification refractory mutation system polymerase chain reaction (ARMS-PCR) was used for genotyping rs11705701, whereas genotyping of rs1470579 and rs2854744 was done employing PCR-restriction fragment length polymorphism (PCR-RFLP) technique. RESULTS: The CC and AA+AC genotypes of rs1470579 conferred an increased risk of PCOS in our population. Regarding the rs2854744, an increased risk of PCOS was observed under the codominant homozygous (TT vs. GG) model by 2.54 fold. The C allele of rs1470579 and T allele of rs2854744 enhanced PCOS risk by 1.97 and 1.46 folds, respectively. Haplotype analysis showed that the Ars1470579Ars11705701 haplotype conferred a decreased risk of PCOS (odds ratio = 0.53, 95% confidence interval = 0.34-0.83, p = 0.006). The AC/GG/GT, AA/GA/GT, AC/GA/GG, and AC/GA/GT genotype combinations of rs1470579/rs11705701/rs2854744 were associated with a decreased risk of the disease. CONCLUSIONS: IGF2BP2 rs1470579 and IGFBP3 rs2854744 enhanced PCOS susceptibility in a Southeastern Iranian population. Further investigation involving larger cohorts representing diverse ethnic backgrounds is needed to confirm the current findings.

A genetic survey of patients with familial idiopathic intracranial hypertension residing in a Middle Eastern village: genetic association study.

BACKGROUND: The aim of this study was to determine whether genetic variants are associated with idiopathic intracranial hypertension (IIH) in a unique village where many of the IIH patients have familial ties, a homogenous population and a high prevalence of consanguinity. Several autosomal recessive disorders are common in this village and its population is considered at a high risk for genetic disorders. METHODS: The samples were genotyped by the Ilumina OmniExpress-24 Kit, and analyzed by the Eagle V2.4 and DASH software package to cluster haplotypes shared between our cohort. Subsequently, we searched for specific haplotypes that were significantly associated with the patient groups. RESULTS: Fourteen patients and 30 controls were included. Samples from 22 female participants (11 patients and 11 controls) were evaluated for haplotype clustering and genome-wide association studies (GWAS). A total of 710,000 single nucleotide polymorphisms (SNPs) were evaluated. Candidate areas positively associated with IIH included genes located on chromosomes 16, 8 (including the CA5A and BANP genes, p < 0.01), and negatively associated with genes located on chromosomes 1 and 6 (including PBX1, LMX1A, ESR1 genes, p < 0.01). CONCLUSIONS: We discovered new loci possibly associated with IIH by employing a GWAS technique to estimate the associations with haplotypes instead of specific SNPs. This method can in all probability be used in cases where there is a limited amount of samples but strong familial connections. Several loci were identified that might be strong candidates for follow-up studies in other well-phenotypes cohorts.

A peek into the population genetics of white teatfish (Holothuria fuscogilva) of Kenya's south coast.

BACKGROUND: The white teatfish, Holothuria fuscogilva, is widely distributed in coastal areas, including waters around coral reefs and seagrasses in the Indo-Pacific. In Kenya, the species is distributed in shallow reefs with higher landings reported from the Vanga-Shimoni-Gazi seascape on the Kenyan south coast. Despite its high exploitation for export and its vulnerable and endangered statuses under IUCN and CITES respectively, Kenya's H. fuscogilva populations and how they may have been impacted by the fishing pressure have not been studied. METHODS: We estimated the genetic diversity and structure of H. fuscogilva population conveniently sampled from three sites in Kenyan south coast using the mitochondrial cytochrome oxidase subunit I (COI) gene sequences. We recorded 30 haplotypes with 43 polymorphic sites across the population. Furthermore, we estimated an overall high haplotype diversity and low nucleotide diversity of estimates of h = 0.970 ± 0.013 and π = 0.010 ± 0.001 respectively. CONCLUSIONS: These preliminary findings suggest several population outcomes, among them a fit population, which require confirming with more comprehensive study to inform strategies for the sustainable exploitation and management of the species.

Biobank-scale inference of multi-individual identity by descent and gene conversion.

We present a method for efficiently identifying clusters of identical-by-descent haplotypes in biobank-scale sequence data. Our multi-individual approach enables much more computationally efficient inference of identity by descent (IBD) than approaches that infer pairwise IBD segments and provides locus-specific IBD clusters rather than IBD segments. Our method's computation time, memory requirements, and output size scale linearly with the number of individuals in the dataset. We also present a method for using multi-individual IBD to detect alleles changed by gene conversion. Application of our methods to the autosomal sequence data for 125,361 White British individuals in the UK Biobank detects more than 9 million converted alleles. This is 2,900 times more alleles changed by gene conversion than were detected in a previous analysis of familial data. We estimate that more than 250,000 sequenced probands and a much larger number of additional genomes from multi-generational family members would be required to find a similar number of alleles changed by gene conversion using a family-based approach. Our IBD clustering method is implemented in the open-source ibd-cluster software package.

HaplotagLR: An efficient and configurable utility for haplotagging long reads.

Understanding the functional effects of sequence variation is crucial in genomics. Individual human genomes contain millions of variants that contribute to phenotypic variability and disease risks at the population level. Because variants rarely act in isolation, we must consider potential interactions of neighboring variants to accurately predict functional effects. We can accomplish this using haplotagging, which matches sequencing reads to their parental haplotypes using alleles observed at known heterozygous variants. However, few published tools for haplotagging exist and these share several technical and usability-related shortcomings that limit applicability, in particular a lack of insight or control over error rates, and lack of key metrics on the underlying sources of haplotagging error. Here we present HaplotagLR: a user-friendly tool that haplotags long sequencing reads based on a multinomial model and existing phased variant lists. HaplotagLR is user-configurable and includes a basic error model to control the empirical FDR in its output. We show that HaplotagLR outperforms the leading haplotagging method in simulated datasets, especially at high levels of specificity, and displays 7% greater sensitivity in haplotagging real data. HaplotagLR advances both the immediate utility of haplotagging and paves the way for further improvements to this important method.

Mitochondrial DNA D-Loop Polymorphisms among the Galla Goats Reveals Multiple Maternal Origins with Implication on the Functional Diversity of the HSP70 Gene.

Despite much attention given to the history of goat evolution in Kenya, information on the origin, demographic history, dispersal route, and genetic diversity of Galla goats remains unclear. Here, we examined the genetic background, diversity, demographic history, and population genetic variation of Galla goats using mtDNA D-loop and HSP70 single-nucleotide polymorphism markers. The results revealed 90 segregating sites and 68 haplotypes in a 600-bp mtDNA D-loop sequence. The overall mean mitochondrial haplotype diversity was 0.993. The haplotype diversities ranged between 0.8939 ± 0.0777 and 1.0000 ± 0.0221 in all populations supporting high genetic diversity. Mitochondrial phylogenetic analysis revealed three Galla goat haplogroups (A, G, and D), supporting multiple maternal ancestries, of which haplogroup A was the most predominant. Analysis of molecular variance (AMOVA) showed considerable variation within populations at 94.39%, evidence of high genetic diversity. Bimodal mismatch distribution patterns were observed while most populations recorded negative results for Tajima and Fu's Fs neutrality tests supporting population expansion. Genetic variation among populations was also confirmed using HSP70 gene fragment sequences, where six polymorphic sites which defined 21 haplotypes were discovered. Analysis of molecular variance revealed a significant FST index value of 0.134 and a high FIS index value of 0.746, an indication of inbreeding. This information will pave the way for conservation strategies and informed breeding to improve Galla or other goat breeds for climate-smart agriculture.

Genotype data for 60 SNP genetic markers associated with eye, hair, skin color, ABO blood group, sex, core Y-chromosome haplogroups in Kazakh population.

OBJECTIVES: The collection of genotype data was conducted as an essential part of a pivotal research project with the goal of examining the genetic variability of skin, hair, and iris color among the Kazakh population. The data has practical application in the field of forensic DNA phenotyping (FDA). Due to the limited size of forensic databases from Central Asia (Kazakhstan), it is practically impossible to obtain an individual identification result based on forensic profiling of short tandem repeats (STRs). However, the pervasive use of the FDA necessitates validation of the currently employed set of genetic markers in a variety of global populations. No such data existed for the Kazakhs. The Phenotype Expert kit (DNA Research Center, LLC, Russia) was used for the first time in this study to collect data. DATA DESCRIPTION: The present study provides genotype data for a total of 60 SNP genetic markers, which were analyzed in a sample of 515 ethnic Kazakhs. The dataset comprises a total of 41 single nucleotide polymorphisms (SNPs) obtained from the HIrisPlex-S panel. Additionally, there are 4 SNPs specifically related to the AB0 gene, 1 marker associated with the AMELX/Y genes, and 14 SNPs corresponding to the primary haplogroups of the Y chromosome. The aforementioned data could prove valuable to researchers with an interest in investigating genetic variability and making predictions about phenotype based on eye color, hair color, skin color, AB0 blood group, gender, and biogeographic origin within the male lineage.

Near-gapless and haplotype-resolved apple genomes provide insights into the genetic basis of rootstock-induced dwarfing.

Dwarfing rootstocks have transformed the production of cultivated apples; however, the genetic basis of rootstock-induced dwarfing remains largely unclear. We have assembled chromosome-level, near-gapless and haplotype-resolved genomes for the popular dwarfing rootstock 'M9', the semi-vigorous rootstock 'MM106' and 'Fuji', one of the most commonly grown apple cultivars. The apple orthologue of auxin response factor 3 (MdARF3) is in the Dw1 region of 'M9', the major locus for rootstock-induced dwarfing. Comparing 'M9' and 'MM106' genomes revealed a 9,723-bp allele-specific long terminal repeat retrotransposon/gypsy insertion, DwTE, located upstream of MdARF3. DwTE is cosegregated with the dwarfing trait in two segregating populations, suggesting its prospective utility in future dwarfing rootstock breeding. In addition, our pipeline discovered mobile mRNAs that may contribute to the development of dwarfed scion architecture. Our research provides valuable genomic resources and applicable methodology, which have the potential to accelerate breeding dwarfing rootstocks for apple and other perennial woody fruit trees.

Positive selection footprints and haplotype distribution in the genome of dromedary camels.

Dromedary camels are a domestic species characterized by various adaptive traits. Limited efforts have been employed toward identifying genetic regions and haplotypes under selection that might be related to such adaptations. These genetic elements are considered valuable sources that should be conserved to maintain the dromedaries' adaptability. Here, we have analyzed whole genome sequences of 40 dromedary camels from different Arabian Peninsula populations to assess their genetic relationship and define regions with signatures of selection. Genetic distinction based on geography was observed, classifying the populations into four groups: (1) North and Central, (2) West, (3) Southwest, and (4) Southeast, with substantial levels of genetic admixture. Using the de-correlated composite of multiple signal approach, which combines four intra-population analyses (Tajima's D index, nucleotide diversity, integrated haplotype score, and number of segregating sites by length), a total of 36 candidate regions harboring 87 genes were identified to be under positive selection. These regions overlapped with 185 haplotype blocks encompassing 1 340 haplotypes, of which 30 (∼2%) were found to be approaching fixation. The defined candidate genes are associated with different biological processes related to the dromedaries' adaptive physiologies, including neurological pathways, musculoskeletal development, fertility, fat distribution, immunity, visual development, and kidney physiology. The results of this study highlight opportunities for further investigations at the whole-genome level to enhance our understanding of the evolutionary pressures shaping the dromedary genome.

Haplotype information of large neuromuscular disease genes provided by linked-read sequencing has a potential to increase diagnostic yield.

Rare or novel missense variants in large genes such as TTN and NEB are frequent in the general population, which hampers the interpretation of putative disease-causing biallelic variants in patients with sporadic neuromuscular disorders. Often, when the first initial genetic analysis is performed, the reconstructed haplotype, i.e. phasing information of the variants is missing. Segregation analysis increases the diagnostic turnaround time and is not always possible if samples from family members are lacking. To overcome this difficulty, we investigated how well the linked-read technology succeeded to phase variants in these large genes, and whether it improved the identification of structural variants. Linked-read sequencing data of nemaline myopathy, distal myopathy, and proximal myopathy patients were analyzed for phasing, single nucleotide variants, and structural variants. Variant phasing was successful in the large muscle genes studied. The longest continuous phase blocks were gained using high-quality DNA samples with long DNA fragments. Homozygosity increased the number of phase blocks, especially in exome sequencing samples lacking intronic variation. In our cohort, linked-read sequencing added more information about the structural variation but did not lead to a molecular genetic diagnosis. The linked-read technology can support the clinical diagnosis of neuromuscular and other genetic disorders.

Phylogeographical and population genetics of Polyspora sweet in China provides insights into its phylogenetic evolution and subtropical dispersal.

BACKGROUND: Geological movements and climatic fluctuations stand as pivotal catalysts driving speciation and phylogenetic evolution. The genus Polyspora Sweet (Theaceae), prominently found across the Malay Archipelagos and Indochina Peninsula in tropical Asia, exhibits its northernmost distribution in China. In this study, we investigated the evolutionary and biogeographical history of the genus Polyspora in China, shedding light on the mechanisms by which these species respond to ancient geological and climatic fluctuations. METHODS: Phylogenetic relationships of 32 representative species of Theaceae were reconstructed based on the chloroplast genome and ribosome 18-26 S rRNA datasets. Species divergence time was estimated using molecular clock and five fossil calibration. The phylogeography and population genetics in 379 individuals from 32 populations of eight species were analyzed using chloroplast gene sequences (trnH-psbA, rpoB-trnC and petN-psbM), revealing the glacial refugia of each species, and exploring the causes of the phylogeographic patterns. RESULTS: We found that Chinese Polyspora species diverged in the middle Miocene, showing a tropical-subtropical divergence order. A total of 52 haplotypes were identified by the combined chloroplast sequences. Chinese Polyspora exhibited a distinct phylogeographical structure, which could be divided into two clades and eight genealogical subdivisions. The divergence between the two clades occurred approximately 20.67 Ma. Analysis of molecular variance revealed that the genetic variation mainly occurred between species (77.91%). At the species level, Polyspora axillaris consists of three lineages, while P. speciosa had two lineages. The major lineages of Chinese Polyspora diverged between 12 and 15 Ma during the middle to late Miocene. The peak period of haplotype differentiation in each species occurred around the transition from the last interglacial to the last glacial period, approximately 6 Ma ago. CONCLUSION: The primary geographical distribution pattern of Chinese Polyspora was established prior to the last glacial maximum, and the population historical dynamics were relatively stable. The geological and climatic turbulence during the Quaternary glacial period had minimal impact on the distribution pattern of the genus. The genus coped with Quaternary climate turbulence by glacial in situ survival in multiple refuges. The Sino-Vietnam border and Nanling corridor might be the genetic mixing center of Polyspora.

Identification and Regulation of Hypoxia-Tolerant and Germination-Related Genes in Rice.

In direct seeding, hypoxia is a major stress faced by rice plants. Therefore, dissecting the response mechanism of rice to hypoxia stress and the molecular regulatory network is critical to the development of hypoxia-tolerant rice varieties and direct seeding of rice. This review summarizes the morphological, physiological, and ecological changes in rice under hypoxia stress, the discovery of hypoxia-tolerant and germination-related genes/QTLs, and the latest research on candidate genes, and explores the linkage of hypoxia tolerance genes and their distribution in indica and japonica rice through population variance analysis and haplotype network analysis. Among the candidate genes, OsMAP1 is a typical gene located on the MAPK cascade reaction for indica-japonica divergence; MHZ6 is involved in both the MAPK signaling and phytohormone transduction pathway. MHZ6 has three major haplotypes and one rare haplotype, with Hap3 being dominated by indica rice varieties, and promotes internode elongation in deep-water rice by activating the SD1 gene. OsAmy3D and Adh1 have similar indica-japonica varietal differentiation, and are mainly present in indica varieties. There are three high-frequency haplotypes of OsTPP7, namely Hap1 (n = 1109), Hap2 (n = 1349), and Hap3 (n = 217); Hap2 is more frequent in japonica, and the genetic background of OsTPP7 was derived from the japonica rice subpopulation. Further artificial selection, natural domestication, and other means to identify more resistance mechanisms of this gene may facilitate future research to breed superior rice cultivars. Finally, this study discusses the application of rice hypoxia-tolerant germplasm in future breeding research.

Male Pedigree Toolbox: A Versatile Software for Y-STR Data Analyses.

Y-chromosomal short tandem repeats (Y-STRs) are widely used in forensic, genealogical, and population genetics. With the recent increase in the number of rapidly mutating (RM) Y-STRs, an unprecedented level of male differentiation can be achieved, widening and improving the applications of Y-STRs in various fields, including forensics. The growing complexity of Y-STR data increases the need for automated data analyses, but dedicated software tools are scarce. To address this, we present the Male Pedigree Toolbox (MPT), a software tool for the automated analysis of Y-STR data in the context of patrilineal genealogical relationships. The MPT can estimate mutation rates and male relative differentiation rates from input Y-STR pedigree data. It can aid in determining ancestral haplotypes within a pedigree and visualize the genetic variation within pedigrees in all branches of family trees. Additionally, it can provide probabilistic classifications using machine learning, helping to establish or prove the structure of the pedigree and the level of relatedness between males, even for closely related individuals with highly similar haplotypes. The tool is flexible and easy to use and can be adjusted to any set of Y-STR markers by modifying the intuitive input file formats. We introduce the MPT software tool v1.0 and make it publicly available with the goal of encouraging and supporting forensic, genealogical, and other geneticists in utilizing the full potential of Y-STRs for both research purposes and practical applications, including criminal casework.

Determination for KIR genotype and allele copy number via real-time quantitative PCR method.

Killer cell immunoglobulin-like receptor (KIR) and human leukocyte antigen (HLA) play crucial roles in regulating NK cell activity. Here, we report a real-time quantitative PCR (qPCR) to genotype all KIR genes and their copy numbers simultaneously. With 18 pairs of locus-specific primers, we identified KIR genes by Ct values and determined KIR copy number using the 2-∆Ct method. Haplotypes were assigned based on KIR gene copy numbers. The real-time qPCR results were consistent with the NGS method, except for one sample with KIR2DL5 discrepancy. qPCR is a multiplex method that can identify KIR copy number, which helps obtain a relatively accurate haplotype structure, facilitating increased KIR research in laboratories where NGS or other high-resolution methods are not available.